Anisotropy of Imbalanced Alfvenic Turbulence in Fast Solar Wind 
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We present the first measurement of the scale-dependent power anisotropy of Elsasser variables in 
imbalanced fast solar wind turbulence. The dominant Elsasser mode is isotropic at lower spacecraft 
frequencies but becomes increasingly anisotropic at higher frequencies. The sub-dominant mode is 
anisotropic throughout. There are two distinct subranges exhibiting different scalings within what is 
normally considered the inertial range. The low Alfven ratio and the different scaling of the Elsasser 
modes suggests an interpretation of the observed discrepancy between the velocity and magnetic 
field scalings, the total energy is dominated by the latter. These results do not appear to be fully 
explained by any of the current theories of incompressible imbalanced MHD turbulence. 



Introduction. The solar wind is an excellent plasma 
turbulence laboratory. Alfvenic fluctuations in the fast 
solar wind (mean velocity > 600 km/s) are well described 
by the incompressible MHD equations despite the colli- 
sionless nature of the plasma The MHD turbulent 
cascade transports energy from large scales to smaller 
scales [e.g. Q until it reaches the ion gyroscale, below 
which another type of turbulence carries energy to yet 
smaller scales [3j-|5[. In the fast wind, the turbulence is 
imbalanced: there is more power in Alfvenic fluctuations 
traveling away from the Sun than toward it @, 0] . There 
is evidence in the slow solar wind Q and from numerical 
simulations Q that balanced turbulence is made up of 
locally imbalanced regions, so understanding imbalanced 
turbulence is probably essential for understanding MHD 
turbulence in general. 

A key property of plasma turbulence is anisotropy 
caused by the magnetic field. Even if the field is not 
strong enough to dominate the thermal pressure, its 
presence makes fluctuations scale differently in the field- 
perpendicular direction than in the field-parallel one, 
with larger power in fluctuations that vary across the 
field. This anisotropy is poorly understood: there is a 
relative dearth of observational data and an abundance 
of mutually contradictory theories. What currently ap- 
pears to be the most compelling theory of the anisotropic 
Alfvenic cascade is based on the assumption of 'critical 
balance' fioj ). to which an assumption of 'dynamic align- 
ment' can be added First posited for balanced cas- 
cades [T3] and later extended to imbalanced ones [l2[ , the 
critical balance conjecture states that the linear and non- 
linear timescales are comparable and predicts anisotropic 
scalings of the fluctuation spectra: E(k±) cx k ± 1 and 



oc /tii . The dynamic alignment conjecture 11 1 
states additionally that the polarizations of magnetic 
and velocity fluctuations align as the energy moves to 
smaller scales, which causes adjustment of the scalings 
to E(k±) cx kj 3 ^ 2 and E(k\\) oc k^ 2 . 

Recent numerical studies have shown a range of seem- 



ingly contradictory behaviors [e.g. 0, which were 

argued to agree, or disagree, with a number of conflict- 
ing theories 0, [T3, [H - 17 1 . So far, observational stud- 



ies of anisotropy have focused on the magnetic field p~8l — 
22j, but a complete analysis of the Alfvenic turbulence 
must include the velocity field. The Elsasser fields [23j ]. 

= V ± B/-y/47rpo, where V and B are the velocity 
and magnetic fields respectively, and po is the average 
mass density, describe oppositely propagating (at local 
Alfven speed) finite-amplitude solutions of the incom- 
pressible MHD equations. There is a good physical case 
for treating these as the primary fields that make up the 
Alfvenic turbulence, which can be thought of as the re- 
sult of their interactions 24\. This Letter presents the 



first in-situ solar wind observation designed to measure 
anisotropy of Elsasser variables as well as the anisotropy 
of the magnetic and velocity fields. The results do not 
appear to be in quantitative agreement with any of the 
current theories. 

Data analysis. We use data obtained by the WIND 
spacecraft at 3-second spin resolution. Magnetic field 
(B) is provided by the MFI instrument, velocity (V) and 
density (p) from the 3DP instrument. Data are taken 
from a long-duration fast stream from days 13 to 20 of 
2008, in which the solar wind speed remained above 550 
km/s for the entire 7-day interval and had an ion plasma 
beta of 1.1. Compressive fluctuations were an order of 
magnitude weaker in power than incompressible, so the 
magnetic fluctuations are dominated by the component 
perpendicular to the mean magnetic field. We denote 
the (dominant) Alfvenic fluctuations traveling away from 
the Sun by Z + for easier comparison to previous work 
[e.g. 0|. Other, shorter fast-stream periods show similar 
results but have larger errors due to insufficient statistics 
in smaller data sets. 
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We use Morlet wavelets [25J to measure the power in 
fluctuations of B and V as a function of time and scale 
El ■ The time resolution of the Morlet wavelet 
is provided by a Gaussian envelope function, the width 
of which changes with scale. This width is used as the 
scale over which we calculate the average density (po) 
and the average angle 9b between the magnetic field and 
the measurement direction (radial) . The power in the El- 
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FIG. 1. Power spectra from WIND data between days 13 and 20 of 2008. (a) Trace of the Fourier and wavelet power spectra 
E v (gray line and blue symbols) and E b (black line and red symbols). The bottom panel shows the Alfven ratio r\ = E v /E b . 
(b) Trace of the Fourier and wavelet power spectra of E + (black line and red symbols) and E~ (gray line and blue symbols). 
The bottom panel shows te = E + /E~ , the high frequency decrease in which may be due to quantization as discussed in the 
text. 



sasser fields is then calculated by combining the wavelet 
coefficients of V and B. In this way we calculate the 
trace power spectra, of B, V, Z + and as a function 
of spacecraft frequency / (proportional to scale under 
Taylor's hypothesis) and the angle 9b- The angle bins 
are 10° wide and the frequency bins are logarithmically 
separated by a factor of 1.6. We stress that anisotropic 
power spectra calculated in this way do not represent 
power in individual components of the vector fields but 
rather the trace power averaged over instances when the 
local mean magnetic field points in the direction given 
by angle 0b with respect to the solar wind velocity. 



The trace power spectra of V and 'B/^/Anpo arc de- 
noted by E v and E b respectively, similarly the power 
spectra of are E^. Subscripts _L and || denote the 
trace spectra corresponding to fluctuations that vary per- 
pendicularly (80° < 9 B < 90°) and parallel (0° < 9 B < 
10°) to the local magnetic field. We refer to the mea- 
sured change in power with ^ at a fixed frequency as 
'power-level anisotropy' and to the change in the mea- 
sured scaling of the power spectrum with 8b as 'spectral- 
index anisotropy'. 

Spectra. We plot E v and E b in the top panel of 
Fig 1(a) The Fourier power of E b is larger than E v 



at all scales, as can be seen in the bottom panel which 
shows the Alfven ratio ta = E v /E b ; at low frequen- 
cies, the magnetic field dominates with ta ~ 0.4, but as 
frequency increases, ta approaches unity. The wavelet 
power measured parallel and perpendicular to the mag- 
netic field is also plotted in Fig 1(a) The magnetic field 



is anisotropic at higher frequencies, with less power in 
the parallel direction. This is similar to the known result 
of magnetic field anisotropy starting at the outer scale 
of turbulence reported for Ulysses data 2(J 21 1, but the 
transition from isotropy to anisotropy occurs at a signif- 
icantly higher frequency than what one would typically 
call the outer scale in such observations (10 -2 Hz rather 
than 1CT 4 Hz). 

The wavelet velocity spectra have lower power than 
E b but behave in a similar manner. As the magnetic 
field becomes anisotropic at frequencies > 10 ~ 2 Hz, so 
does the velocity, with the parallel power decreasing and 
the perpendicular power dominating. Thus, E b and E v 
track each other, suggesting a common source. Error bars 
are not plotted here to keep the figures clear; errors on 
individual wavelet power observations grow as frequency 
decreases, being the same size as the circle and square 
markers at the highest frequency and increasing in size to 
about 50% of the wavelet value at the lowest frequencies 
(~ 10~ 3 Hz). 



Fig 1(b) shows the Fourier power spectra E + and E 



(top panel) calculated for the same period as Fig 1(a) 
The Elsasser variables show imbalance, with E + dom- 
inant over E~ at all frequencies. This is quantified in 
the bottom panel which shows te = E + /E~ , at all fre- 
quencies for this period ve > 0; for purely anti-sunward 
Alfven waves we expect te = oo. Fig 1(b) also shows the 



has equal power in both directions for / < 10 Hz but 



Elsasser power spectra parallel and perpendicular to the 
local mean magnetic field. The dominant (Z + ) modes 
are isotropic in power at frequencies lower than 10~ 2 
Hz, but grow increasingly anisotropic at higher frequen- 
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FIG. 2. Power and spectral index anisotropy of the two subranges of the power spectra in Fig. 1. (a) Subrange 1 (10 _1 and 
10 -2 Hz): variation of power (top panel) at 3.26 x 10 -2 Hz and spectral index (bottom panel) with 8b- (b) Subrange 2 (10 -2 
and 10~ 3 Hz): variation of power (top panel) at 3.1 x 10 -3 Hz and spectral index (bottom panel) with Ob- 



cies. This behavior is similar to that for the magnetic 
field in Fig 1(a) which is expected since ta < 1- The 



weaker (Z J mode behaves differently, with power-level 
anisotropy measured at all scales. 

Thus, we have found that what is usually thought of 
as the 'inertial range' of solar wind turbulence in fact 
consists of two distinct subranges, each about a decade 
wide in this data interval: the higher-frequency Subrange 

1 (~ 1CP 2 - 10 _1 Hz) and the lower-frequency Subrange 

2 (~ 1CP 3 - 10~ 2 Hz). The scalings and anisotropy are 
clearly different in these two subranges. We now proceed 
to quantify these differences in terms of power-level and 
spectral-index anisotropy. 

Scaling and anisotropy in Subrange 1. The top panel 
of Fig 2(a) shows the power at / = 3.26 x 10~ 2 Hz, 



plotted against 9b for each of the four fields. The bot- 
tom panel shows the spectral indices measured over Sub- 
range 1. There is power- level anisotropy in all 4 variables 
with the ratios E b ± /E b = 2.4 ± 0.2, E\/E\\ = 2.2 ± 0.2, 



E+/E+ = 2.4 ± 0.2, 



and E'/ET; = 2.0 ± 0.2, although 



these ratios change with scale for E b , E v and E + due to 
their spectral-index anisotropy. Errors in power are the 
standard deviation of the data contributing to each mean, 
and errors in the spectral index are calculated from the 
least-squares fit of a straight line to the wavelet power 
spectra on a log-log scale. 



The bottom panel of Fig 2(a) attests that the spectral 
index of E + is almost identical to E b , with steeper spec- 
tra at smaller 9 b and the spectral index tending to —2 at 
9b = 0, consistent with standard theories based on criti- 



cal balance 1(J For larger angles (8b > 15°), E + and 
E b have spectral indices around —5/3. The E v spectral 
indices are slightly shallower than the E b and E + values 



throughout and tend towards —3/2 for 0b > 15°, consis- 
tent with the shallower gradient of the E v Fourier spec- 



trum compared to E (see Fig 1(a) and [e.g. |26j). E is 



clearly different from the other fields: it shows no obvious 
spectral-index anisotropy which also has much shallower 
values around —1.3. Thus, the sub-dominant Elsasscr 
field has power-level anisotropy but its anisotropy does 
not change with scale. 

There is a caveat concerning the E~ results, however. 
The WIND 3DP velocity data are subject to quantiza- 
tion; the velocity observations are digitized in such a way 
that the velocity appears measured in discrete steps at 
small scales rather than as a smoothly varying time se- 
ries. A careful analysis shows that the effect on E v and 
its anisotropy is small due to the low power of the quanti- 
zation; however the Elsasser variables are more strongly 
affected. The quantization acts to decorrelate the mag- 
netic and velocity fluctuations at frequencies higher than 
/ ~ 3 x 10~ 2 Hz, increasing the power in the weaker 
Elsasser variable with a corresponding slight decrease 
in the dominant mode. Only the effect on the weaker 
Elsasser mode is significant relative to the errors here, 
and as such the results for / > 3 x 10 -2 Hz represent 
a minimum bound on the power-level anisotropy of the 
weaker Elsasser variable. The quantization is unlikely 



to have affected the power- level anisotropy in Fig 2 (a 



since we have chosen a scale where the noise effect is 
small. The E~~ spectral indices are shallower than they 
should be over this frequency range and the quantization 
may remove any spectral-index anisotropy. The power- 
level anisotropy does not change significantly at lower fre- 
quencies, however, confirming the lack of spectral-index 
anisotropy in general, if not the exact values of the spec- 
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tral indices in this subrange. 

Scaling and anisotropy in Subrange 2. The top panel 
of Fig |2(b)| shows the power of each variable plotted vs. 
6 B at / = 3.1 X 1(T 3 Hz. E+, E b and E v are ap- 
proximately flat, exhibiting no systematic power-level 
anisotropy. E~ , however, retains a level of anisotropy 
similar to that seen in Subrange 1, with greater power at 
larger angles than in the parallel direction. The power- 
level anisotropy ratios are E\/E^ = 1.1 ± 0.3, E v ± /E^ = 
0.7 ±0.3, E+/E+ = 0.9 ±0.3, and E];/E^ = 1.7 ±0.4. 
The bottom panel shows the spectral indices plotted 
against 8b- Subject to larger errors than in subrange 1, 
the spectral indices of E b and E + do not show any mea- 
surable change with angle, and hover around the —5/3 
and —3/2 values. E v shows similar spectral indices to 
E b and E + for 8 b < 45° but its spectral index gets shal- 
lower at larger angles, increasing to around —1.3 in the 
perpendicular direction. The spectral index of E~ shows 
an entirely different trend: it is flatter at around —1.5 
in the parallel direction, becoming steeper at larger an- 
gles until it is less than —2 at 90°. This is in agreement 
with the result that the sub-dominant Elsasser mode has 
a steeper spectrum than the dominant one at lower fre- 
quencies [a • 

Conclusions and discussion. We have presented 
the first observations of power-level and spectral-index 
anisotropy of Elsasser variables in imbalanced fast solar 
wind turbulence. The two remarkable conclusions of this 
study are as follows, (i) What is usually thought of as a 
simple 'inertial range' is in fact split into two subranges. 
In the data used here the split is at / ^ 10~ 2 Hz, where 
both the scaling behavior and the anisotropy of all fields 
change, (ii) While V, B and Z + have similar anisotropy 
(very little in the lower-frequency subrange, increasing 
with frequency in the higher-frequency subrange), the 
sub-dominant Elsasser field is completely different: it has 
power-level anisotropy but in a scale-independent way at 
higher frequencies, and scales steeply in the perpendic- 
ular direction at lower frequencies. Z + and B behave 
very similarly at all scales and the anisotropy of B and 
V appear related; the change from isotropy to anisotropy 
occurs at the same frequency in both, although the spec- 
tral indices are shallower for the velocity than for the 
magnetic field. 

We believe that our results shed some light on the 
theoretically puzzling observations that suggest different 
spectral scalings of V and B, —3/2 for the former and 
—5/3 for the latter J26|, a behavior never observed in sim- 
ulations or envisioned in theories. In the solar wind, the 
Alfvcn ratio is significantly below 1 resulting in B and Z + 
scaling in the same way. Since Z~ has a shallower scal- 
ing, the V scaling lies between these two behaviors and 
the closeness of the spectral exponent to —3/2 might, in 
fact, be a coincidence. The dominant contribution to the 
energy comes from Z + (and so B), which has a robust 
—5/3 scaling. 

Although these results in their entirety are not pre- 



cisely reproduced by any of the extant simulations or 
theory, there are some points of similarity. The MHD 
simulations of [IH show different scaling and anisotropy 
in the two Elsasser variables for imbalanced turbulence 
(the dominant Z + scales in a critically balanced way and 
the weaker Z~ scales differently, although the scaling of 
Z~ does not agree with our findings). The theory of 
[Til l also predicts different scaling for Elsasser variables, 
with the concept of 'pinning' of E + and E~ to the same 
power at the dissipation scale. In particular, the results 
from a mixture of weak and strong turbulence in (l6| are 
qualitatively similar to the results here, with a break in 
the middle of what appears to be an inertial range and a 
change from steeper to shallower scaling of Z , although 
again the anisotropy of Z does not agree with our re- 
sults. 
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